Elliptic flow and energy loss of heavy quarks in ultra-relativistic heavy ion collisions 
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The space-time propagation of heavy quarks in ultra-relativistic heavy ion collisions is studied 
within the partonic transport model Boltzmann Approach of MultiParton Scatterings (BAMPS). 
In this model heavy quarks interact with the partonic medium via binary scatterings. The cross 
sections for these interactions are calculated with leading order perturbative QCD, but feature a 
more precise Debye screening derived within the hard thermal loop approximation and obey the 
running of the coupling. Within this framework the elliptic flow and the nuclear modification factor 
of heavy quarks are computed for RHIC and LHC energies and compared to available experimental 
data. It is found that binary scatterings alone cannot reproduce the data and, therefore, radiative 
corrections have to be taken into account. 

PACS numbers: 25.75.-q, 25.75.Bh, 25.75.Cj, 12.38.Mh, 24.10.Lx 



I. INTRODUCTION 

In ultra-relativistic heavy ion collisions a medium is 
produced that behaves like a nearly perfect fluid PH3, 
i.e. has a small viscosity to entropy density ratio. This 
medium is thought to be a new state of matter, which 
consists of quarks and gluons and, therefore, is called 
quark-gluon plasma (QGP). 

Heavy quarks, in particular, are an ideal probe for this 
medium. Due to their large mass, they are produced in 
the early stage of the collision, where the energy den- 
sity is large [5[. Consequently, they traverse the QGP 
and interact with the rest of the medium. Due to these 
processes, their distributions arc modified and can reveal 
- via experimentally accessible observables like the el- 
liptic flow and nuclear modification factor - information 
about the properties of the medium. Furthermore, heavy 
quarks are rare and tagged by their flavor, which renders 
them as an unique probe, even after hadronization, due 
to flavor conservation. 

The experimentally measured elliptic flow vi and nu- 
clear modification factor Raa of heavy flavor electrons 
are comparable to that of light hadrons. This re- 
sult was surprising since it was thought that the radiative 
part of the energy loss of heavy quarks is suppressed due 
to the "dead cone effect" @, [lOj]. Whether the large el- 
liptic flow and strong suppression is due to collisional or 
radiative interactions - or both (or even other effects) - 
is under investigation 

04H. 

In the present study we investigate the contribution of 
clastic scatterings with a running coupling and an im- 
proved Debye screening inspired by hard thermal loop 
calculations to «2 and Raa- 

This article is organized as follows. In the next sec- 
tion we present our model BAMPS. Section ITTT1 explains 



the modifications we employed to the standard leading 
order cross section. In Sec. IIVI we show our results for 
RHIC and LHC and compare them to the experimental 
data where possible. Finally, we conclude with a short 
summary. 



II. PARTON CASCADE BAMPS 

The Boltzmann Approach of MultiParton Scatterings 
(BAMPS) [13, HH is a 3 + 1-dimensional partonic trans- 
port model, which solves the Boltzmann equation, 
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for on-shell partons and perturbative QCD (pQCD) in- 
teractions. It includes elastic and also inelastic gluonic 
(g) interactions, gg gg, ggg gg, and gg -!• ggg, the 
last one being important, for instance, for thermalization 
[HHlj], elliptic flow [Hill], or jet-quenching [HHl of 
gluons. 

Heavy quarks (Q) are produced in initial hard parton 
scatterings or in the QGP in the reaction gg — > QQ (the 
back reaction of this process is also possible but negligible 
at RHIC and LHC). Studies on heavy quark production 
within the QGP at RHIC and LHC energies have been 
carried out with BAMPS and can be found in Refs. 0, 
[27l - [29j . In the present model heavy quarks interact with 
the gluonic medium via elastic collisions, gQ — > gQ or 
gQ — > gQ. The implementation of radiative corrections 
such as gQ — > gQg is planned for the future. 

The incorporation of light quarks (q) is underway and 
will be presented soon. However, the impact on heavy 
quarks will be very small, since gluons dominate the early 
stages of the heavy ion collision and the cross section of 
qQ — > qQ is suppressed compared to gQ — > gQ due to a 
smaller color factor. 

The initial heavy quark distributions are generated 
with the Monte Carlo event generator for next-to-leading 
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Figure 1: Differential invariant cross section of heavy flavor 
electrons as a function of transverse momentum for proton- 
proton collisions with ^J~s = 200 GeV simulated with mc@nlo 
without shadowing effects. For comparison experimental data 
8, 35] arc also shown. 

order (NLO) calculations MC@NLO [3(1 HJ. To com- 
pare the initial distributions to the experimental data 
from proton-proton collisions, the heavy quarks are frag- 
mented to D and B mesons, which, consequently, decay 
into heavy flavor electrons. 

The Peterson fragmentation function [32J of a heavy 
quark to a heavy meson H 

N 

D h/q(z) = — -2 (2) 



is used for the fragmentation process with N being 
a normalization constant, z — \ph\I\vq\, and cq — 
0.05 (0.005) for charm (bottom) quarks. The decay of 
heavy flavor mesons into electrons is performed with 
PYTHIA 8.1 [HI Hj|. Figure [T] compares measured heavy 
flavor electrons to our initial distributions obtained with 
mc@nlo for a fact orization a nd renormalization scale of 
HP = (x R = 0.65 + Ml for charm (M c = 1.3 GeV) 
and up = = 0.4 \Jp\ + M 2 for bottom quarks 
(Mb = 4.6 GeV). These (in principle, arbitrary) scales 
are chosen in such a way that a good agreement with the 
experimental data is found. 

As a note, large theoretical uncertainties in the heavy 
quark distributions exist due to uncertainties in the par- 
ton distribution functions, renormalization and factoriza- 
tion scale, and heavy quark masses |5|, 36]. In particular, 
the relative contributions of charm and bottom quarks to 
the electron spectrum are not fully settled yet. Although 
we choose the scales such that the electron spectrum at 
RHIC is reproduced and employ the same parameter set 
for LHC initial conditions, there are considerable uncer- 
tainties in the initial heavy quark distributions at LHC. 

The initial distributions for nucleus-nucleus collisions 
are obtained by scaling from proton-proton collision with 



the Glauber model. In the present study nuclear effects 
for the parton distribution functions, for instance, shad- 
owing or the Cronin effect, are not taken into account for 
heavy quark production since the impact on the heavy 
quark distributions at intermediate and high transverse 
momentum pt is rather small (ll| . 

The initial gluon distributions are obtained either with 
pythia or the mini-jet model. More details can be found 
in Refs. [HH. 

Since charm and bottom quarks are very rare probes 
(there are about four heavy quark pairs at mid-rapidity 
in a central Au+Au collision at RHIC compared to about 
800 gluons Q), one needs to simulate several million 
events to yield sufficient statistics to obtain their spec- 
tra and elliptic flow. However, the most time consuming 
part of these event simulations is the computation of the 
interactions among the sizeable number of gluons, for 
which we do not need as many statistics. Therefore, we 
generate a smaller sample of pure gluonic events and use 
them several times as a background medium in which 
heavy quarks generated with mc@nlo are placed. This 
treatment is in perfect accordance with the conventional 
BAMPS model. The only difference to full simulations 
is the neglect of medium modifications induced by the 
heavy quarks, which is a very good approximation given 
the small number of heavy quarks in a heavy ion collision. 



III. MODIFICATION OF THE CROSS SECTION 

The cross sections for Q + g — > Q + g can be calculated 
in leading-order (LO) pQCD. To treat this accurately, we 
explicitly take the running of the coupling into account. 
Furthermore, we employ a more precise Debye screening 
compared to the standard procedure. This is done by 
comparing the energy loss calculated with the LO cross 
section and the hard thermal loop (HTL) approximation 

MM- 

The LO cross section [38[ diverges for small Man- 
delstam t due to the gluon propagator in the t chan- 
nel. In thermal field theory, however, long-range in- 
teractions (which correspond to small t) are screened 
by the medium, an effect that originates formally from 
thermal loop corrections. The resummation and renor- 
malization of these loop corrections results in a screen- 
ing of the gluon propagator with its thermal self-energy 
ntherm(w, q) and replacing the bare coupling with the 
running coupling a s (t) [33, [H| 



t t - n t horm(w,<?) 

with (w,q) = Pi — P% (the four- momentum difference 
of the incoming and outgoing heavy quark in the center- 
of-mass frame). Nevertheless, calculations with the self- 
energy are too involved to be incorporated in a transport 
model. However, we approximate the self-energy by a 
screening mass /i 2 = nrn^, which is proportional to the 
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The Debye mass is calculated by [2( 
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where N c = 3 denotes the number of colors, v g — 16 the 
gluon degeneracy and / the distributions of gluons and 
quarks. In equilibrium and with Boltzmann statistics this 
simplifies to m 2 D cq = ^-(N c + nj)T 2 . For consistency 
we employ a running a s = a s (t) also in the calculation 
of the Debye mass. 

The prefactor k in Eq. |(4)| is mostly set to 1 in the 
literature without a sophisticated reason. However, one 
can fix this factor analytically by comparing the energy 
loss per unit length dE/dx of the LO cross section in- 
cluding n (see Eq. |(4)[ ) to the energy loss within the hard 
thermal loop approach [H, |37[ • 

In the following we will outline the calculation for a 
constant coupling a s since, in this case, k can be deter- 
mined analytically. The generalization to a running cou- 
pling is given in Ref. [17] . The heavy quark collisional 
energy loss within the HTL approximation was calculated 
in Refs. [l7], HO, IH| for quantum statistics. Analogously, 
for Boltzmann statistics 1 the collisional energy loss of a 
heavy quark with energy E and mass M in a thermalized 
medium with temperature T is given in the high energy 
limit by 



dE 

da; 

with 
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h ~ 0.251 n f + 0.747, 



(6) 



(7) 



where g = ln2-l/4-7/3 and ft = (31/9) ln2-101/108- 
II7/9 with 7 being the Euler-Mascheroni constant and 
rif the number of flavors. 

The energy loss calculation above is done within the 
HTL approach. An analogous calculation can be re- 
peated with the screened (according to Eq. |(4)[ ) LO cross 
section instead of the full HTL calculation. One ends up 
with a result which is very similar to Eq. |(6)[ 
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Only the argument of the logarithm in the first term 
changes. From this, one can read off, that the screening 
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Figure 2: Energy loss per unit length dE/dx of a heavy quark 
jet in a static and thermalized medium of gluons (n / = 0) with 
temperature T = 300 MeV as a function of the jet energy E. 
The numerical result from BAMPS with the LO cross section, 
a constant coupling a s = 0.3, and a screening mass with 
K — 0.184 is compared to the analytical result from Eq. (6) 
for Boltzmann statistics. In addition, the curve for quantum 
statistics [13, Bll] is also shown. 



prefactor k has to be k = l/(2e) f=s 0.184 ps 0.2 to obtain 
the same result as in the HTL approximation [3 l37j . 
The determination of k for a running coupling yields also 
a result close to this value (T^ |. 

In Fig. [5] the numerical result with BAMPS and k = 
0.184 is compared to the analytic formula [(6)[ For large 
jet energies, that is, the regime in which the analytic for- 
mula is valid, the agreement with the numerical result is 
very good. For comparison the analytic curve for quan- 
tum statistics is shown as well, which is about 25 % larger 
than the Boltzmann curve for E ~ 10 GeV. 

The running coupling is evaluated at the renormaliza- 
tion scale While complete results are /in independent 
due to renormalization group flow equations, approxi- 
mate results do have a residual /in dependence. The re- 
sulting uncertainty can be reduced if one chooses a scale 
relevant for the physical problem at hand. In the t chan- 
nel, for instance, large logarithms of t occur in next-to- 
leading order due to vacuum contributions to the self- 
energy and vertex corrections. These logarithms can be 
absorbed in a s via renormalization by choosing p 2 R = t 
[i~7L l39l | . An analogous line of argument holds also for 
the s and u channels using their characteristic momenta 
instead of t [13, H3]- Thus, we evaluate a s (fi 2 R ) at the 
virtuality of the respective channel, that is, s — M 2 , t, 
and u — M 2 for the s, t, and u channels, respectively. 

An effective description of the running coupling can be 
obtained from measurements of e + e~ annihilation and 
non-strange hadronic decays of r leptons and continued 
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to the time-like region [1J, |37J, |42 



a s (Q 2 ) 



Air L 



Q 2 < 



(9) 

\ - 7r~ 1 arctan(L + / 7 r) Q 2 > w 



with /3 = 11 - f n/ and L± = ln(±Q 2 /A 2 ) with A = 
200 MeV. If a s (Q 2 ) is larger than a' s nax = 1.0 it is set 
to a™ ax . This cutoff procedure is in line with a soft 
average of a s given by the universality hypothesis [37], 
|43j]. We checked that our results of the energy loss are 
not very sensitive on the exact value of A or a™ ax since 
the contribution from the soft part of the process is very 
small. 
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IV. ELLIPTIC FLOW AND NUCLEAR 
MODIFICATION FACTOR 

The elliptic flow 



1>2 




(10) 



(p x and p y are the momenta in the x and y directions 
with respect to the reaction plane) and the nuclear mod- 
ification factor 



Raa 



d 2 N AA /dp T dy 
N hhl d 2 N pp /d PT dy 



(11) 



are suitable observables to study the energy loss of heavy 
quarks and their coupling to the medium. The larger the 
i>2 the stronger heavy quarks interact with the medium 
and adopt its momentum anisotropy. A small value of 
Raa indicates a strong suppression and, therefore, large 
energy loss of heavy quarks in the medium. Unfortu- 
nately, D and B mesons that stem from charm and 
bottom quarks, respectively, cannot be reconstructed di- 
rectly at RHIC yet. Therefore, one measures single elec- 
trons that are decay products of D and B mesons. The v-i 
and Raa of these heavy flavor electrons are comparable 
to the respective values for light hadrons 

To compare to the experimental data we also perform 
the fragmentation of heavy quarks to D and B mesons 
and the consecutive decay to electrons. The final distri- 
bution of the electrons is very similar to that of heavy 
quarks, only shifted to lower transverse momentum pt- 
For small pr Peterson fragmentation is not a good de- 
scription of the hadronization process and other schemes 
like coalescence might be more appropriate. However, 
we postpone the implementation of more sophisticated 
hadronization schemes in BAMPS to the future. 

For elastic heavy quark scatterings with the medium, 
the LO pQCD cross section with a constant coupling and 
a standard Debye screening is too small to explain the 
measured data j27|. However, the improvements pre- 
sented in Sec. IIIII increase the cross section and yield 
results which are much closer to the experimental data. 
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Figure 3: (a) Elliptic flow V2 and (b) nuclear modification 
factor Raa of heavy quarks and heavy flavor electrons with 
pseudo-rapidity [77] < 0.35 for Au+Au collisions at RHIC with 
an impact parameter of b = 8.2 fm. The curves are obtained 
with pythia and mini-jet initial conditions (IC) for the glu- 
ons. The cross section of gQ — > gQ is multiplied with the 
factor K = 4. To estimate the uncertainty of this K factor 
we plotted for mini-jet IC the electron curves for K = 3 and 
K = 5 as gray bands. For comparison, data of heavy flavor 
electrons for the centrality class of 20 %— 40 % Q are shown. 



In Fig.[3]we show the elliptic flow and nuclear modifica- 
tion factor of heavy quarks and heavy flavor electrons for 
a non-central Au+Au collision at RHIC obtained with 
BAMPS at the end of the QGP phase, that is, after 
the energy density of a given cell has dropped below 
0.6GeV/fm 3 and interactions are not allowed any more 
23]. For these curves the cross section of gQ — > gQ has 
been multiplied with an artificial factor K = 4 to be com- 
patible with the data. BAMPS studies on the energy loss 
of light partons reveal that radiative contributions are 
dominant for light particles [25l |26| . For heavy quarks 
these contributions are expected to be suppressed due to 
the dead cone effect, but it has been shown that they are 
still slightly larger than the elastic contributions [lj, [l5[ . 
Consequently, we assume that the implementation of ra- 
diative corrections like gQ — > gQg could account for the 
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missing factor of 4. However, it needs to and will be 
checked in a forthcoming study whether these contribu- 
tions have indeed the same effect as a constant K factor. 
The implementation of NLO processes for heavy quarks 
will complement 2f>3 interactions for gluons, which are 
already present in BAMPS [201 ] . Furthermore, the con- 
sideration of quantum statistics would also enhance the 
cross section as is shown in Fig. [5] for the energy loss and 
can, therefore, explain part of the missing factor of 4. 

Of course, a constant factor K — 4 is only an esti- 
mate for radiative contributions and quantum statistics. 
To explore the uncertainty of this factor the curves for 
K = 3 and 5 are also plotted in Fig. [3] as gray bands 
for electrons and mini-jet initial conditions. Both ob- 
servables t>2 and Raa are not very sensitive on the exact 
value of K, although K = 4 gives the best agreement 
with the data. 

Figure [3] compares the elliptic flow and heavy quark 
suppression obtained with two different initial conditions 
for the gluons. Both the pythia and the mini-jet sce- 
nario lead to comparable curves, especially for larger px- 
This indicates that both heavy flavor observables are not 
very sensitive on the initial light parton distributions as 
long as the elliptic flow of light hadrons is reproduced, 
which has been shown for mini-jet initial conditions in 
Refs. [22|-[24j]. In addition, we confirmed this also with 
pythia initial conditions and will present the results in 
a forthcoming paper. However, in Ref. [44| it was found 
that differences in the medium evolution can lead to mod- 
ifications of the heavy quark suppression and elliptic flow 
up to a factor of 2. 

With K = 4 the agreement with the data is very good 
for large px- We emphasize that both V2 and Raa are 
described simultaneously within the same partonic trans- 
port model. For small px Peterson fragmentation is not 
suitable and coalescence might play an important role, 
which modifies the vi and Raa due to the contribution 
of light quarks. 

These results are obtained with initial heavy quark 
distributions from mc@nlo. In addition, we employed 
heavy quark initial conditions from NLO calculations of 
Refs. [H, HI] and found a good agreement with the curves 
shown in Fi g, p In contrast to pythia which we used 
in Refs. [H, l29j] for the initial heavy quark distributions, 
the electron spectrum from MC@NLO is slightly steeper 
for low pt which results in a less suppressed Raa and, 
therefore, better agreement with the experimental data 
for pt < 3 GeV. 

In previous studies 0, [27], HsJ we have found that at 
LHC energies a sizeable fraction of the produced charm 
quarks is created during the evolution of the QGP and 
not only in the initial hard parton scatterings. These 
secondarily produced charm quarks can, of course, con- 
tribute to the V2 and Raa of heavy quarks. However, 
the impact on these observables is very small and only 
affects the region of low px since secondary charm quarks 
are produced with small momenta. Due to their early 
production time which is usually less than lfm/c [5[ 



secondary charm quarks have enough time to interact 
with the medium, lose energy, and build up elliptic flow. 
Consequently, also in the low px region, V2 is barely 
changed by secondary charm quarks. The fraction of 
bottom quarks produced in the medium is so small that 
the assumption that all of them are created in the initial 
hard parton scatterings is justified [5j. In short, since 
the contribution of secondarily produced heavy quarks is 
insignificant we neglect it in the following. 

In addition to the uncertainties of the initial heavy 
quark distributions (see Sec. |H]), there are also sizeable 
uncertainties concerning the bulk medium at LHC. For 
initial gluon conditions from pythia we obtain a final 
dNS/dy fa 1770 and dE T /dy 1570 GeV of gluons at 
mid-rapidity in central collisions after the energy density 
has dropped below 0.6GeV/fm 3 for every cell. Unfor- 
tunately, experimental data for dEx/dy is not available 
yet and it is rather involved to compare the measured 
diV c h/d77 of charged hadrons to the number of gluons 
due to non-perturbative effects at hadronization. How- 
ever, we checked that the ratio of diV c h/d?7 at LHC [46| 
and dN c h/dr] at RHIC [47| for central collisions, (without 
taking errors into account) 
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diVch/dT?) 



LHC 



dN ch /d V \ 



RHIC 



1601 
687 



2.33 



(12) 



is nearly equal to the ratio of the final dN 9 /dy at LHC 
and RHIC obtained with BAMPS, which is for pythia 
initial conditions 



R ._ dN3/dy\ LHC 
9 ' dN9/dy\ Rmc 



1770 
740 



2.39 



(13) 



This indicates that the scaling to LHC has been per- 
formed accurately. Nevertheless, a comparison to mea- 
sured dEx/dy will be an additional test. 

For mini-jet initial conditions one can vary the mo- 
mentum cutoff po for Pb+Pb collisions at the LHC to 
obtain an R g which is close to R c h- This procedure leads 
to po = 3.5 GeV and 



Ra 



1820 
800 



2.28 . 



(14) 



The final transverse energy distribution for this cutoff is 
dE T /dy « 1810 GeV. 

For an impact parameter of b = 8.2 fm the final val- 
ues of gluons at mid-rapidity are dN 9 /dy s» 450(440) 
and dE T /dy « 410 GeV (570 GeV) for pythia (mini-jet 
with po — 3.5 GeV) initial conditions. The ratio of the 
former to our RHIC value of dN^/dy w 210(200) for 
b = 8.2 fm and pythia (mini-jet) initial conditions is 
R g = 2.14 (2.20) which is close to the experimental ratio 
of LHC and RHIC of i? c h = 2.37 for the centrality class 
20%-30% m S3. 

Figure 2] shows and Raa of initially produced heavy 
quarks for a Pb+Pb collision at LHC with ^/snn = 
2.76 TeV and an impact parameter of b = 8.2 fm. As 
for RHIC (cf. Fig. \S§ we employ a cross section for the 
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Figure 4: As in Fig.O but for Pb+Pb collisions at LHC with y/sJ^N = 2.76 TeV. 



heavy quark interactions with gluons, which incorporates 
the running coupling, the improved Debye screening, and 
K = 4. Both for heavy quarks and heavy flavor electrons 
the elliptic flow for PYTHIA and mini-jet initial conditions 
agree very well. Heavy quarks and thus also heavy flavor 
electrons are slightly less suppressed for mini-jet initial 
conditions compared to pythia initial conditions. This 
effect can also be seen for gluons and is probably due to a 
smaller initial gluon number at mid-rapidity for the mini- 
jet case compared to the PYTHIA scenario, which results 
in a less opaque medium before the QGP is chemically 
equilibrated. 

The elliptic flow and suppression of electrons from 
charm quarks are slightly larger than at RHIC. However, 
the contribution from bottom quarks becomes important 
also at smaller px- As a consequence, the vi and Raa at 
LHC are very similar to the RHIC results. 

V. CONCLUSIONS 

We presented the BAMPS results of the elliptic flow 
and nuclear modification factor of heavy flavor electrons 
at RHIC and LHC. For this study, elastic interactions 
with the partonic medium are taken into account in LO 



pQCD. The cross section of these processes obey the 
running of the coupling and include a Debye screening 
motivated from HTL calculations, which is more precise 
compared to previous approaches. However, elastic scat- 
tering alone does not lead to the sizeable elliptic flow 
and nuclear suppression measured at RHIC. To yield the 
same values as the data, the cross section must be mul- 
tiplied with a factor K = 4. We assume that radiative 
contributions to the interaction can account for this phe- 
nomenological factor, which we will check in an upcoming 
study. The elliptic flow and nuclear modification factor 
at LHC are found to be of the same order as at RHIC. 
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